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Abstract 

We investigate a new model for the finite one-dimensional quantum oscillator based upon 
the Lie superalgebra st(2|l). In this setting, it is natural to present the position and momen- 
*-£h I turn operators of the oscillator as odd elements of the Lie superalgebra. The model involves a 

parameter p (0 < p < 1) and an integer representation label j. In the (2j + l)-dimensional repre- 
[ sentations Wj of sl(2|l), the Hamiltonian has the usual equidistant spectrum. The spectrum of 

the position operator is discrete and turns out to be of the form iVfc, where k = 0, 1, . . . ,j. We 
construct the discrete position wave functions, which are given in terms of certain Krawtchouk 
polynomials. These wave functions have appealing properties, as can already be seen from their 
plots. The model is sufficiently simple, in the sense that the corresponding discrete Fourier 
transform (relating position wave functions to momentum wave functions) can be constructed 
! explicitly. 

m 

1 Introduction 

Finite quantum harmonic oscillator models (simply referred to as finite oscillator models) are of 
importance in optical image processing [7J, or in models where only a finite number of eigenmodes 
can exist such as in signal processing [2rS]- Quantum kinematics of finite oscillators has also been 
used to remove infinities or divergences in quantum theory or quantum field theory [22]. The main 
^ . idea underlying these finite oscillator models is to replace the Heisenberg algebra of the standard 
quantum oscillator, which allows only infinite-dimensional representations, by a "closely related" 
algebra with the same dynamics but which allows finite-dimensional representations. 

More concretely, for a one-dimensional finite oscillator one considers three (essentially self- 
adjoint) operators: a position operator q, its corresponding momentum operator p and a Hamilto- 
nian H which is the generator of time evolution. These operators should satisfy the Hamilton-Lie 
equations (or the compatibility of Hamilton's equations with the Heisenberg equations): 

[H,q] = -ip, [H,p]=iq, (1) 

in units with mass and frequency both equal to 1, and h = 1. The third relation of the canonical 
oscillator, [q,p] = i, is dropped since otherwise the corresponding algebra (the Heisenberg algebra) 
has infinite-dimensional representations only. Instead, one requires [5]: 

• all operators q, p, H belong to some (Lie) algebra (or superalgebra) A; 

• the spectrum of H in (unitary) representations of A is equidistant. 
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The most popular model is with A = su(2) (or its enveloping algebra), see [SHI]- This model is also 
underlying the theory in [22j. The relevant finite-dimensional representations are the common su(2) 
representations labelled by an integer or half-integer j. In such a representation, the Hamiltonian is 
taken as H = Jq+j + \, where Jo = J z is the diagonal su(2) operator, and thus the spectrum of H is 
n + ^ (n = 0, 1, . . . , 2j). The operators q and p are linear combinations of the other su(2) operators 
J + and J_ , such that the relations (Q} are satisfied. They turn out to have a finite spectrum given by 
{—j, — j + • • • , +j} [5j, see Figure 1(a). In this context, one can also construct the discrete position 
and momentum wave functions. For the 5u(2) case, these are given by Krawtchouk functions 
(normalized symmetric Krawtchouk polynomials), i.e. Krawtchouk polynomials with parameter 
p = 1/2. These discrete wave functions have many interesting properties, and their shape is 
reminiscent of those of the canonical oscillator [5J . 

In |13y i4j. deformations of su(2) were considered by introducing a deformation parameter a > 
— 1. For the even-dimensional representations |13| (j half- integer) , this led to an alternative model 
of the finite oscillator with the spectrum of H again given by n + \ (n = 0, 1, . . . , 2j), and with the 
spectrum of the position operator q given by 

±(a + k) (k = l,2,...,j + -), 

i.e. a finite and mostly equidistant spectrum apart from a gap of size 2a + 2 in the middle, see 
Figure 1(b). For the odd-dimensional representations [13] (j integer), this led to a similar model 
with the same spectrum of H but with the spectrum of q given by (Figure 1(c)) 

0, ±y/k(2a + k + l), (k = l,...,j). 

In both deformations, the position wave functions could be constructed explicitly and turned out 
to be related to normalized Halm (or dual Hahn) polynomials [13[I14]. It was also shown how these 
wave functions could be interpreted as the finite version of a paraboson oscillator. 

Recall however that Wigner [25] was among the first to drop the canonical commutation relation 
[q,p] = i, proposing a model with relations (pQ) with the extra condition 

H=\{f+q 2 ). (2) 

This is known as the Wigner quantum oscillator (or paraboson oscillator) (19|.I20| . This paraboson 
oscillator still has an infinite energy spectrum of the form n + a (with a positive representation 
parameter a), and the spectrum of the position operator is the real axis. The canonical oscillator 
is recovered from the paraboson oscillator in the representation with a = 1/2, i.e. one of the 
representations of the paraboson oscillator coincides with the canonical oscillator jl2 |, llT |. [T8] . The 
algebraic structure equivalent with (pQ) and (|2|) is the Lie superalgebra osp(l|2) [T21[TS]. In this 
context, the position and momentum operators q and p are odd elements of the Lie superalgebra, 
whereas H is an even element. This observation, and the fact that the canonical oscillator fits in 
one of the osp(l|2) representations, leads to the idea that it is perhaps more natural to consider 
the position and momentum operators of alternative oscillator models as odd elements of a Lie 
superalgebra rather than just (even) elements of a Lie algebra. 

Following this last idea, we propose and investigate in this paper a finite oscillator model based 
on the Lie superalgebra sl(2|l). Indeed, apart from osp(l|2), sl(2|l) is the simplest Lie superalgebra 
that can be considered as a superversion of the Lie algebra su(2) [2l] . The Lie superalgebra st(2| 1) 
has moreover a class of representations Wj of dimension 2j + 1 (j integer) which are similar to the 
common su(2) representations, see Section 2. 

In the finite oscillator model studied here in Section 3, the Hamiltonian H is a diagonal operator 
with spectrum n + \ (n = 0, 1, . . . , 2j) in Wj. The position operator q is an arbitrary (self-adjoint) 
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odd element from s[(2|l), and the form of the momentum operator p follows from (pQ); the model 
introduces in a natural way a real parameter p with < p < 1. Our investigation shows that the 
spectrum of q in the representation Wj is still very simple, given by 

±Vk (k = Q,l,...J), 

see Figure 1(d). In order to prove this statement, we construct the eigenvectors of q in Wj explicitly. 
For this construction, we need some properties of Krawtchouk polynomials (but now with general 
parameter p), a well known set of discrete orthogonal polynomials. 

The properties of the position and momentum wave functions for the new sl(2|l) oscillator 
model are investigated in Section 4. In particular, we study first some plots of the wave functions 
for different values of the parameter p. These discrete wave function plots are rather different from 
the 5u(2) plots in [5]; only for small p- values these plots show some similarity with the plots of the 
continuous canonical oscillator wave functions. We explore also the behaviour of the discrete wave 
functions when the representation parameter j tends to infinity, and discover a relation with the 
paraboson oscillator. 

When both position wave functions and momentum wave functions are known explicitly, one 
can examine the transformation that relates the two. In the canonical case this is just the Fourier 
transform. In a finite oscillator model, this gives a discrete version of the Fourier transform, 
determined by a Fourier matrix F. For the sl(2|l) model, we manage to give an explicit expression 
for an arbitrary matrix element F^i of F, see Section 5. The properties of the matrix F are studied 
and prove to be similar to those of the standard Discrete Fourier Transform used in spectral analysis. 
Our paper closes with some concluding remarks in Section 6. 

2 The Lie superalgebra s((2|l) and a class of representations 

The Lie superalgebra gl(2|l) is well known: it can be defined as the set of all 3 x 3 matrices x = (xij) 
with rows and columns labelled by indices i,j = 1,2,3. As a basis in g[(2|l) we choose the Weyl 
matrices e^, i, j = 1, 2, 3, where the odd elements are {en, en\i = 2, 3}, and the remaining elements 
are even. The Lie superalgebra bracket is determined by 

[ eij , e kl j = eij e kl - (-l)deg(^)deg(e fc! ) ewe . . = ^ _ deg(e„)^ efc . (3) 

Note that the even part g[(2|l)o = 0t(2) ©fll(l), where gl(2) = span{ey|i, j = 1,2} and fll(l) = 
span{e33}. For elements x of g[(2|l), one defines the supertrace as str(x) = x\\ + X22 — %33- The 
Lie superalgebra g[(2|l) is not simple, and one can define the simple superalgebra sl(2|l) as the 
subalgebra consisting of elements with supertrace 0. As a basis for sl(2|l), it is convenient to follow 
the choice of |10j . where one can find the actual matrices of the basis [10\ p. 261]: 

F + = e 32 , G + = ei 3 , F~ = e 31 , G~ = e 23 , (4) 

1 1 
H = -(en - e 22 ), E + = e 12 , E~ = e 2 i, Z = -(e n + e 22 ) + e 33 . (5) 

So s[(2|l) has four odd (or 'fermionic') basis elements F + , F~ , G + , G~ and four even (or 'bosonic') 
basis elements H, E + , E~ , Z . The basis for the gl(2) subalgebra is {H, E + , E~} and the gl(l) = 
U(l) subalgebra is spanned by Z. The basic Lie superalgebra brackets can be found in [TOj, p. 261] 
or, in a different notation, in [21j . For the odd elements, the anti-commutators are given by 

{F ± ,G ± } = E ± , {F ± ,GT} = Z T H; 

{F ± ,F ± } = {G ± ,G ± } = {F ± ,F^} = {G ± ,G^} = 0. (6) 
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For the even elements, the commutators are given by 

[H,E ± ] = ±E ± , [E+,E~]=2H, [Z, H] = [Z, E±] = 0. (7) 
The mixed commutation relations read: 

[H,F ± ] = ±±F ± , [Z,F ± ] = ±F ± , [£±^±1 = 0, [E*,F±] = -F*; 

[H,G ± ] = ±^G ± , [Z,G±]=-±G±, [^,^1=0, [^, G ±] =G T. ( 8 ) 

The irreducible representations of sl(2|l) have been studied by Scheunert et al [21] and Marcu [16] ; 
for a summary, see [TOj §2.53]. The superalgebra sf(2| 1) has typical and atypical irreducible rep- 
resentations. Here, we shall consider a class of atypical irreducible representations, labelled by a 
non-negative integer j (these are denoted by iv-(j/2) in [10]). In order to describe these represen- 
tations explicitly, let us choose a basis for the representation space Wj of the form 

\j,m), m = -j, -j + 1,. . . (9) 

So dim(Wj) = 2j + 1. For the actions of the st(2| 1) basis elements on these vectors, it is handy to 
use the following "even" and "odd" functions, defined on integers n: 

£(n) = 1 if n is even and otherwise, 

0(n) = 1 if n is odd and otherwise. (10) 

Note that 0(n) = 1 — £(n), but it is convenient to use both notations. The actions of the odd 
generators are now given by: 

F ± \j,m) = ±0(j-m)^ J± ™ + 1 \j,m±l), 

G±\j,m) = ±£(j - m)^J J -^ \j,m±l). (11) 

The actions of the even generators can in principle be computed from (|6j), and are 

Z\j,m) = -£(j -m) 3 - \j,m) - 0(j - m)(^— ) \j,m), 
m 

E ± \j,m) = -£(J - m)y / (j=Fm)(j±m + 2) \j,m± 2) 

+ ~0(j - mW(j T m - l)(j ± m + 1) \j, m ± 2). (12) 

It is easy to verify that the relations ©-([HI) are satisfied for these actions. One can also see that with 
respect to the even subalgebra g[(2)©gl(l) = su(2)ffif7(l), Wj decomposes as (|; — 1)©(^-; — ^-), 
where (I; b) denotes the su(2) © U(l) representation "with isospin I and hypercharge b" [TO] . 

The above representation is a star representation (or unitary representation) for the adjoint 
operation 

tf = Z, H^ = H, (^ ± )t = ^ ) (F±)t = _ G =F j ( G ±)t = _ F T ( 13) 
compatible with the positive definite inner product on the representation space Wf. 

(j,m\j,m') = 8 m ,m'- ( 14 ) 
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3 The sl(2|l) model for a finite one-dimensional oscillator 



We wish to investigate how sl(2|l) can be used for a finite oscillator model. Inspired by the seminal 
paper [7J on the su(2) model for a finite oscillator, and by the requirements for the spectrum of the 
Hamiltonian H, one should take 

H = 2H + j + ~ (15) 

as operator for H in the representation space Wj. This operator is diagonal, self-adjoint, and has 
indeed the equidistant spectrum: n + ^ (n = 0, 1, 2, . . . , 2j). Next comes the choice for the position 
operator q, which should be a self-adjoint element of sl(2|l). There are two reasons to choose an 
odd element of st(2|l). First, choosing an even element would yield a model that is essentially the 
same as the su(2) model of [TJ. Second, as explained in the introduction, in Wigner's paraboson 
oscillator the position and momentum operator are elements of the odd part of the Lie superalgebra 
osp(l|2) |12 y i8 1 l25]. And since the canonical oscillator is a special case of the paraboson oscillator 
(corresponding to one particular osp(l|2) representation |12|. I18j ) . the position and momentum 
operator can also be considered as odd elements of a superalgebra in a particular representation. 
The most general real self-adjoint odd element of st(2|l) is given by 



AF + + BG + - BF~ - AG- 



(16) 



with A and B real constants. An overall constant does not play a crucial role, so let us assume 
that there is some normalization like A 2 + B 2 = 1. Consider now first the case that A and B have 
the same sign, say both positive (the case A positive and B negative will be very similar, and is 
described at the end of this section). In that case, one can write 



q = ^F + + ^T^G + - ^T^F--^G-, (0<p<l). 

We shall in fact consider < p < 1 and later view the values p = and p = 1 as a limit. 
Once q is fixed, the form of p follows from the first equation of (pQ), and thus 



(17) 



P 



i{y/p F + + y/T^p G + + yjl^p F~ + ^pG~ 



(0<p< 1). 



(18) 



With these operators, ([I]) and the conditions described in Section 1 are satisfied, and we can truly 
speak of an sl(2|l) model for the oscillator. 

Now it remains to study the operators q and p, in particular their spectrum and eigenvectors 
in the representation Wj. Note that, due to the actions (fTT]) . in the (ordered) basis {\j,j}, \j,j — 
1), . . . , |j, —j + 1), \ j, —j)} of Wj, the operators q and p are tridiagonal matrices. In particular, for 
q, one has: 



where 
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The matrix form M p of p is similar. For these matrices, we need to study the spectrum and the 
eigenvectors. It is at this point that Krawtchouk polynomials play a role. Krawtchouk polynomials 
K n (x;p, N) of degree n in the variable x, with parameter p are defined by [HUHIT5]: 

K n (x;p,N) = 2 F 1 (~ n \ T X ;-), (n = 0, 1, . . . , N) (21) 



-N p 

in terms of the hyper geometric series 2 F± [H[8j[23] (which is terminating here because of the ap- 
pearance of — n as numerator parameter). Their (discrete) orthogonality relation reads [Hlll|IT5]: 

N 

Y^w{x;P,N)K l {x;p,N)K n (x;p,N) = h{n;p,N)5i n , (0 < p < 1) (22) 

x=0 

where 

N\^ f , sN—x ,_ »rv ,/_._ ,n n!(iV-n)! fl-pV 



u,(x;p,iV) = ( Ip^l-p)""* (* = 0,l,...,iV); h(n;p,N) = " v ^ 7 ' ( — ) .(23) 
For the orthonormal Krawtchouk functions we use the notation: 



K n (x;p,N) ee ^™l X '' P >y K n (z-,p,N). (24) 
y/i(n;p,iV) 

Now we are able to describe the eigenvalues and (orthonormal) eigenvectors of M q . In this context, 
two sets of Krawtchouk polynomials play a role: K n (x;p,j) (with N = j) and K n (x;p,j — 1) (with 
iV = j-l). 

Proposition 1 Let M q (i.e. the matrix representation of q) be the tridiagonal (2j + 1) x (2j + 1)- 
matrix (|19p and let U = (JJkl)o<k,K2j be the (2j + 1) x (2j + 1) -matrix with matrix elements: 

U 2nJ = (-l) n K (n;p,j), ne{0,l,...,j}; U 2n+1;j = 0, n G {0, ... ,j — 1}; (25) 
(-l) n - 

U^nJ-k = U 2n ,j+k = -f= K k (n;p,j), n £ {0, 1, . . . , j}, fc 6 {1, . . . , j}; (26) 

(-l) ra ~ 

U 2n+1J - k = -U 2n+1 j +k = -f=-Kk-i(ri]p,j - 1), n G {0, 1, . . . , j - 1}, fe € {1, . . . , j}. (27) 

T/ien C7 is an orthogonal matrix: 

UU T = U T U = I. (28) 
The columns ofU are the eigenvectors of M q , i.e. 

M g U = UD, (29) 

where D is a diagonal matrix containing the eigenvalues of M q : 

D = diag(- y 7 !, - VJ 3 !, ■ • • , -V2, -1, 0, 1, y/2, . . . , vT 3 !, y/j). (30) 



Proof. Using the orthogonality of the Krawtchouk polynomials, and the explicit expressions (|25|) - 
(|27|) . a simple computation shows that (U t U)m = Ski- Thus U T U = I, the identity matrix, and 
hence UU T = I holds as well. 
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It remains to verify (|29|) and that the eigenvalues are indeed ([30]) , Due to the tridiagonal 



form (|T9|) of M q , one has: 

[M q Uy 2n k = S n U2n-l,k + Rn+\U2n+l,k, (31) 
(M q U) 2n+l k = R n +\U2n,k + £n+1^2n+2,fc- (32) 

For the first case (f3Tj) . we need to consider three distinct subcases, according to k belonging to 
{0,1,..., j- 1}, to {j + l,j + 2,..., 2j} or k = j. For k G {0, 1, . . . ,j - 1}, this gives: 



(M q U)2n,j-k — S n U2n~l,j-k + -Rn+1^2n+l 



-1)" n-^.^. ^ , ("I) 



n+1 



y/1 -j>yjnk k _ x {n - l;p,j — 1) H 7= — y^W nk k _ x {n;p,j - 1) 



72 v rv "" 1V ' V2 

__ (_i)n+i ( v ^)"+ 1 (yT^?)J-"- 1 v ^T)! 

\/2 \/ n! 0' - ^)!/i(A; - j - 1) 

1 — p 

x [{j ~ n)K k _i(n;p,j - 1) - n( )iffc-i(n - l;p,j - 1)]. 

P 

For the last linear combination between squared brackets, the backward shift operator formula for 
Krawtchouk polynomials [T5J (9.11.8)] can be applied, and yields 



v2 \/nl{j - n)lh{k - l;p,j - 1) 
= -VkU 2n , J -k=(UD) 2nj _ k . 

For the other two subcases with first index In, the computations are similar. For the case (|32p . we 
need again to consider three subcases. Now one finds for k € {0, 1, . . . ,j — 1}: 

{M q U)2n+l,j-k = Rn+lU2n,j-k + S n +lU2 n +2,j -k 



(-If 

(-!)«+! ( v ^)" +1 (v / r Tr py" n vj! 



\n+l 



V^VJ - nK k (n;p,j) H -7= — yl -p\/n + l^(n + l;p,j) 

[K k (n + l;p,j) - K k (n;p,j)]. 



y/2 \/nJ(J -n - l)\h(k;p,j) 

For the last linear combination, the forward shift operator formula for Krawtchouk polynomials |15] 
(9.11.6)] can be applied, and yields 

{MqU) 2 n+l,j-k = 7^ . ... =F= , -^-l «iP,J ~ 1) 

V2 y/n\(j -n-l)\h(k;p,j)P3 
= ~Vk U 2n +i lj - k = {UD) 2n+l j _ k . 

The other two subcases are similar. This completes the proof. □ 
The above proposition gives, besides the eigenvectors, also the spectrum of the position operator 
q in the representation Wj. It will be appropriate to denote these g-eigenvalues by q k , where 
k = -j, -j + l,...,+j, so 

q± k = ±y/k, k = 0, 1, . . . , j. (33) 

To our knowledge, this is the first time we come across a tridiagonal operator with such a spectrum. 
Note also that the spectrum of q is independent of the parameter p in ()1T[) . but the eigenvectors 
themselves do depend on p. 
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As far as the eigenvectors of q are concerned, these are the columns of the matrix U. It will be 
useful to introduce a notation for these eigenvectors: the orthonormal eigenvector of the position 
operator q in Wj for the eigenvalue q k , denoted by \j, is given in terms of the standard basis by 



\j,Qk)= ^2 U j+™,j+k\j,-m). 



(34) 



m=-j 



For the matrix representation M p of p, the analysis is essentially the same. Without going into 
the details, we give the final result: 

Proposition 2 Let M p be the tridiagonal {2j + 1) x (2j + \)-matrix representing p in Wj, and let 
V = (Ykl)o<k l<2j be the (2j + 1) x (2j + 1) -matrix with matrix elements 



V 2 u = -i{-l) k U 2kth 



V- 



2fc+lJ 



-D k u, 



2k+l,U 



(35) 



where U is the matrix determined by (|25p - (|27p . Then V is a unitary matrix, VV^ = V^V = I. 
The columns ofV are the eigenvectors of M p , i.e. 



M„V = VD, 



(36) 



where D is the same diagonal matrix as in Proposition 1. In other words, the eigenvalues of p are 
also given by: 

- V7, - x/T 7 !, • • • , -V% -1,0, l,y/2,..., v^T 7 !, sfj- (37) 



The matrix V of eigenvectors satisfies: 



V T V 



( o 


V -i 



-1 \ 
-1 

/ 



and 



V = JU where J = -i diag(i°, i 1 , i 2 , i 3 , . . . , i 2j ) = diag(-i, 1, i, -1, . . .)• 



(38) 



(39) 



The last assertions follow from the explicit expressions (|35p . (|25p - (|27p and the orthogonality prop- 
erties of Krawtchouk polynomials. 

Also here, it will be useful to denote the ^-eigenvalues by p k , where k = —j, — j + 1, . . . , +j (so 
p±k = ±\/fc, k = 0, 1, . . . , j), and to write the normalized eigenvectors as: 



\j,Pk) = Yl V 3+m,j+k\j,- 



(40) 



m=-j 



We end this section with two remarks. First, let us briefly return to the remaining case (|16p 
with A positive and B negative; or, without losing generality, A = ^Jp and B = —y/l — p. The 
corresponding matrix M' q is then the same as in (jl9>(P0"|) . but with S k replaced by -^/T^pVk. 
In other words, one can write M' q = D\M q D\ where D\ = diag(l, 1, —1, —1, 1, 1, —1, —1, 1, 1, . . .). 
This implies that M' q U' = U'D, where D is the same matrix as in Proposition 1 and U' = D\U . To 
conclude for this second case: the eigenvalues remain the same, and the matrix U' of eigenvectors 
is the same as that of the first case, up to sign changes in rows. For this reason, we shall not return 
to this second case, and just continue with the first case (I17h for our analysis. 

Secondly, we have so far considered (fT7|) with < p < 1, but what about the cases p = and 
p = 1? In these limiting cases, the form of M q in (I19p remains valid, and since the eigenvalues 
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of M q in (|30p are independent of p, the eigenvalues are again given by (|30p. For the matrix of 
eigenvectors U, one can simply take the right limit p — > or the left limit p — )■ 1 in the matrix [/ 
given in Proposition 1. For example, it is easy to compute: 



/ 



U = lim U, 

p>0 



1 

V2 









V2 











-1 





1 








1 





1 





-1 











1 


1 











1 



\ 



V 



(41) 



/ 



where the general form of the orthogonal matrix Uo (with only two nonzero elements per row, 
starting from the second row onwards) is clear from the above. The left limit p — > 1 yields a similar 
matrix form for U. 



4 Position and momentum wave functions and their properties 

The position (resp. momentum) wave functions of the sl(2|l) finite oscillator are the overlaps 
between the g-eigenvectors (resp. ^-eigenvectors) and the //-eigenvectors (or equivalently, the Jo- 
eigenvectors \j,m)). Let us denote them by ^fj m (q) (resp. tp^_ m (p) ), where m = j, j — 1,. .. , —j, 
and where q (resp. p) assumes one of the discrete values qk (resp. pk) (k = —j, + 
Observe that we have denoted the momentum variable of the wave function by p in order not to 
confuse with the parameter p (0 < p < 1) which appears in expression (117j) and which is also the 
parameter of the Krawtchouk polynomials occurring here. In the notation of the previous section 
(and where we want to emphasize the dependence on the parameter p), we have 

vflrrMk) = (j, q k ) = U j+mJ+k , (42) 

^fl m {Pk) = {j, -m\j,p k ) = V j+m j +k . (43) 

Let us consider the explicit form of these wave functions, first for the position variable. For j + m 
even, j + m = 2n, and for positive g-values one has 

^(q k ) = i^K k (n;p,j), n = 0,l,...,j, k = l,...J, (44) 
or equivalently, with q k = y/k {k = 1, 2, . . . , j): 



The expression for <p^(q_ k ), with q_ k = —\J~k~ {k = 1, 2, . . . , j) is also given by the right hand side 
of (|45p . So 4>2n * s an even function. For the argument 0, one simply has 



When j + m is odd, j + m = 2n + 1, one finds for positive q- values 
V2 



< \ — lr ~ 

^2n+i(Qk)= y -^=-K k ^(n;p,j-l), n = 0, 1, . . . ,j - 1, k=l,...,j, (47) 
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or, with q k = y/k (k = 1,2, ... 



4ixfe) = -^0 - l)y n|(j _ x _ n)]{k _ _ k)] 2^1 [ _ J + 1 ,-)■ (48) 

The expression for (ft^n+iiQ-k), with = — \/k (k = 1, 2, . . . ,j) is given by minus the right hand 

side of (USD, and <^£ +1 (0) = 0. So 4n+i is an odd function. 

Let us now consider some plots of these discrete wave functions. In Figure 2 we have plotted 
these functions for the representation j = 10 (so discrete plots with 2j + 1 = 21 points). We have 
considered three values for the parameter p: p = 0.1, p = 1/2 and p = 0.9. The spectrum of q is 
independent of p, so for each of these three cases we plot points corresponding to the values ±\/k 
(k = 0, 1, . . . , 10) on the horizontal axis. For each of the considered p-values, we have plotted (p^ 
(the ground state), (ffi (the first excited state), 4>2 and ^3^- 

The behaviour of these discrete wave functions is reminiscent of that of the corresponding 
continuous wave functions of the canonical oscillator, especially when p is small. As p increases, 
the wave function values for positions near the origin tend to decrease. In fact, for increasing 
p-values, the behaviour of the discrete wave functions rather tends to the corresponding wave 
functions of the paraboson oscillator (see e.g. Figure 3 of [13]). 

Secondly, it is interesting to investigate what happens when the representation parameter j 
increases, i.e. when the dimension of the representation Wj increases. For this purpose, we have 
plotted the ground state and the first excited state, for a fixed p-value, and for the values j = 10, 
j = 30 and j = 60 in Figure 3. These plots remind of the shape of (continuous) paraboson wave 
functions ^n(q) for increasing values of a. In Figure 4 we have illustrated some of these functions, 
and the similarity is indeed striking. This can also be confirmed by a limit calculation. The limit 
is not simply 

hm 4f)( % ); 

that would just yield zero, as the non-zero contributions are shifted further away from the origin 
(see Figure 3). Instead, we need to involve another set of orthogonal polynomials, the dual Hahn 
polynomials R n (X(x);'j,5, N) defined by [El (9.6.1)] 

R n (X(x); 1 ,5,N)= sF 2 (~ n ~ X ^ + ^ + Y + 1 -,l), (n = 0,1, . . . ,N), (49) 

where A(x) = x(x + 7 + <5 + l). These polynomials satisfy a discrete orthogonality relation (7 > — 1, 
<5>-l): 

TV 

w(x; 7, 5, N)R m (X(x);j, 5, N)R n (X(x); 7, 5, N) = h(n; 7, S, N)5 mn , (50) 

x=0 

with w{x; 7, 5, N) and h{n; 7, 5, N) given by [151 (9.6.2)]. Let us also fix a notation for the orthonor- 
mal functions: 

RnMx); 7) 6, N) = J ^hMl ^x^, 5, N). (51) 
y h(n; 7, 5, N) 

Consider now some positive parameter a > 0. Let j be the representation parameter, n = 0,1,... ,j, 
and < p < 1. The following limit is obvious: 

>im ^f^-*L* + *; + 1 ;l)=^f-^:* i i). (52) 
a^oo V -j^2pa + l ) \ -j pi 
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(53) 
(54) 



In other words, 

lim Rn(\(k);2pa,2(l -p)a,j) = K n (k;p,j). 

a— >oc 

Using some elementary limits for the expressions in (I51jl . one can in fact show: 

lim Rn(X(k);2pa,2(l -p)a,j) =K n (k;p,j). 

a— >oo 

For the wave functions (|45p under consideration, this means 

( \ (— l) n ~ (—l) n ~ i—X\ n 

</>%>(q k ) = y —jl-K k (n;p,j) = ^-K n (k;p,j) = ±-jL- Jim i^(A(fc); 2pa, 2(1 -p)a,j), (55) 

for values q k = V~k, k = l,2,...,j. On the other hand, one finds with [13^ Eq. (30)] that for 
jx 2 = A(Jfe) = k(k + 2a + 1), 



/-n, -A:, + 2a + 1 
hm 3 F 2 [ ; 1 

j^oo \ —j, 2pa + 1 

-(a) 



1 2pa + l' 



{2pa + l) n 



4 2 ^)(x 2 ), (56) 



where L n is a (generalized) Laguerre polynomial. Following the limit computations in |13|. §4], 
one obtains 

(— l) n 

lim ^Lj 1 ^R n {ix 2 ;2pa,2{l-p)a,j) 



(-iy 



IV. 



T(n + 2pa + 1) 



(57) 



where ^^\x) is the paraboson wave function with parameter a (see |13} (A. 11)]). Combining (I55D 
and ([57]) . it follows indeed that for large j-values the behaviour of the wave functions ^nG?*) i s 
the same as the behaviour of the paraboson wave functions ^^rf" ^( x ) f° r large values of a. 

For the explanation above, we have used even wave functions; clearly for wave functions of 
degree 2n + 1 the computation is similar and the conclusion is the same. 

5 The corresponding discrete Fourier transform 

In canonical quantum mechanics, the momentum wave function (in L 2 (R)) is given by the Fourier 
transform of the position wave function (and vice versa): 



ip{p) 



1 



e- ipq 4>{q)dq. 



In the present situation we are dealing with discrete wave functions, and an analogue of this should 
be viewed as follows. Let 



(P) ( 



Qk) 



( $\<tk) \ 



( 4\pk) \ 



\ V#(p*) J 



(k = -j,...,+j). 



(58) 



In this case, the corresponding discrete Fourier transform is defined as the matrix F = (i<fc/)_,-<fc z<+j 
relating these two wave functions: 



k=-j 



(59) 
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By (|42p -(j43p. the columns of U consist of the column vectors <f>^(pk) (k = —j--, +j) and similarly 
for the matrix V . So (|59p means that V = U ■ F, or: 

F = U T -V = U T JU, (60) 

with J given by (]39h . Using the explicit matrix elements for U and J , this leads to the following 



form of the matrix elements of F 

i 
2 



F i-k,jTi = F j+kJ±i = ~^S(k,l;p,j) ±~S(k-l,l- l;p,j - 1), (k, I = 1, . . . , j); (61) 



V2 



S(k,0;p,j), (k = l,...,j); 



Fjj = — iS(0, 0;p, j), 



where 



S(k,l;p,j) = ^(-l) n K fe (n;p,i)^(n;p, j). 



(62) 
(63) 

(64) 



n=0 



This last expression is easy to simplify, using [JH Proposition 3] or Eq. (12), p. 85]. One finds: 

-k,-l 

-j 



S(k,l;p,j) 




P)) 



(k+Q/2, 



i - 2 P y 



-k-l 



2*1 



1 



(65) 



4p(l — p) 

So from (|65p and ()61|) -(|63 [) we have an explicit form for the elements F^i of the matrix F. Just as 
in [H], one can state the following properties of the discrete Fourier transform matrix F: 

Proposition 3 The (2j + l) x (2j + 1) -matrix F is symmetric, F T = F , and unitary, F^F = FF^ = 
1. Furthermore, it satisfies F 4 = I, so its eigenvalues are ±1, ±i. A set of orthonormal eigenvectors 
of F is given by the rows of U , determined in Proposition d The multiplicity of the eigenvalues 
depends on the parity of j. When j = 2n is even, then the multiplicity of — i, 1, i, —1 is n + 1, n, n, n 
respectively. When j = 2n + 1 is odd, then the multiplicity of — i, 1, i, — 1 is n + 1, n + 1, n + 1, n 
respectively. 

Proof. The symmetry of F is easily seen from the expressions (16ip - (l63p . The unitarity of F 
follows from (|60p . the orthogonality of the real matrix U and J^J = I. Again using ()60|) and the 
orthogonality of U, one finds F 2 = FF T = V T UU T V = V T V. But the explicit form of V T V is 
known, see (I38p . Since (V T V) 2 = I, the result F A = I follows. So the eigenvalues can only be 
±l,±i. From the last part of (|60p . one has FU T = U T J . So the columns of U T (or the rows of 
U) form a set of orthonormal eigenvectors of F, and the eigenvalues of F are found in the diagonal 
matrix J , see ([39]) . □ 

We are dealing here with a simple discrete Fourier transform matrix F, with parameter < p < 
1: the matrix elements of F are given by 2F1 expressions (|65p . the eigenvectors of F are given by 
the rows of U in Proposition 1, and the eigenvalues of F (with their multiplicities) are given above. 
F still has the property that it transforms position wave functions into momentum wave functions. 

Note that due to expression (i65j) . the matrix F is particularly simple when p = 1/2, as most of 
its matrix elements are zero. For j = 3, the form of F for p = 1/2 reads 



2F 



and it is clear how this generalizes for arbitrary j. 



( 








1 


-iV2 


-1 








\ 







1 


— i 





— i 


-1 









1 


— i 











— i 


-1 






-iy/2 

















-\y/2 






-1 


— i 











— i 


1 









-1 


— i 





— i 


1 







V 








-1 


-iy/2 


1 








/ 



(66) 
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6 Conclusions 



We have, for the first time, explored the possibility of using a Lie superalgebra as the basic structure 
underlying a finite oscillator model. This was inspired by the idea that the position and momentum 
operators of an oscillator model could most naturally be represented by odd (rather than even) 
elements of a superalgebra. In the case presented here, we have taken the Lie superalgebra st(2| 1) 
for this purpose, being a simple generalization of the Lie algebra su(2). 

The s[(2|l) representations most suitable to use in the model are the (2j + l)-dimensional 
representations Wj. Indeed, in the standard basis for these representations the Hamiltonian is 
diagonal and the position and momentum operators are self-adjoint tridiagonal matrices. The 
most general form for the position operator q involves a parameter p (0 < p < 1), see (I17p . The 
first main result of the paper is the determination of the eigenvalues and eigenvectors of q in 
explicit form (Proposition 1). The spectrum is very simple, see ()33|) : the eigenvectors are in terms 
of Krawtchouk polynomials K n (x;p,j) or K n (x;p,j — 1). The eigenvalues and eigenvectors of the 
momentum operator p are similar, see Proposition 2. 

The matrix U of eigenvectors of q is interesting from a second point of view. Indeed, its rows 
correspond to the discrete position wave functions of the finite oscillator. These wave functions 
have been examined in Section 4, by means of plots and by investigating a limit. Although the 
discrete wave functions are also given in terms of Krawtchouk polynomials, as in the su(2) case, 
the behaviour is quite different. This is because in the su(2) case [5j the Krawtchouk polynomials 
appearing in the wave functions are K n (x;^,2j + 1), whereas here we have a combination of 
K n (x;p,j) and K n (x;p,j — 1). For large values of j, the discrete wave functions tend to certain 
paraboson wave functions. 

Since the position and momentum wave functions have fairly simple forms, we are able to 
construct the matrix F that transforms them into each other. This matrix is a discrete analogue of 
the Fourier transform. The matrix elements of F are terminating %Fi series. The discrete version 
of the Fourier transform F has many properties in common with the standard Discrete Fourier 
Transform, see Proposition 3. 

Apart from the finite-dimensional representations Wj considered here, the Lie superalgebra 
sl(2|l) also has an interesting class of infinite-dimensional representations in ^ 2 (Z + ). It would be 
tempting to study sl(2|l) oscillator models in these representations, and we hope to tackle this 
problem in a future paper. 
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((') ••••••• • ••••••• 

(el ) ••••••• • ••••••• 

Figure 1: Plots of a typical spectrum of the position operator: (a) for the su(2) model [5]; (b) for 
the u(2) a model [33]; (c) for the su(2) a model [H]; (d) for the sl(2|l) model. 
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Figure 2: Plots of the discrete wave functions 4>i (<?) in the representation with j = 10, for 
n = 0, 1, 2, 3 and for p = 0.1 (left column), p = 0.5 (middle column) and p = 0.9 (right column). 
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Figure 4: Plots of the continuous paraboson wave functions (x) for n = (left column) and 
n = 1 (right column), and with a varying: a = 5, 10, 30. 
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